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ABSTRACT 

Observations of ultraluminous infrared galaxies (ULIRGs) with an achieved 
resolution approaching the diffraction limit in the mid-infrared from 8-25 //m 
using the Keck Telescopes are reported. We find extremely compact structures, 
with spatial scales of < 0.3" (diameter) in six of the seven ULIRGs observed. 
These compact sources emit between 30% and 100% of the mid-infrared energy 
from these galaxies. 

We have utilized the compact mid-infrared structures as a diagnostic of 
whether an AGN or a compact (100 - 300 pc) starburst is the primary power 
source in these ULIRGs. In Markarian 231, the upper limit on the diameter of the 
12.5 //m source , 0.13", shows that the size of the infrared source must increase 
with increasing wavelength, consistent with AGN models. In IRAS 05189-2524 
and IRAS 08572-1-3915 there is strong evidence that the source size increases with 
increasing wavelength. This suggests heating by a central source rather than an 
extended luminosity source, consistent with the optical classification as an AGN. 

The compact mid-infrared sources seen in the other galaxies cannot be used 
to distinguish the ultimate luminosity source. If these ULIRGs are powered 
by compact starbursts, the star formation rates seen in the central few hundred 
parsecs far exceed the global rates seen in nearby starburst galaxies, and approach 
the surface brightness of individual clusters in nearby starburst galaxies. 

Subject headings: luminous infrared galaxies, infrared, galaxies individual: 
Arp 220; Markarian 231; Markarian 273; UGC 5101; IRAS 05189-2524; 
IRAS 08572+3915; IRAS 17208-0014 
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1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs; defined as tfiose systems witfi luminosity 
from 8 - 1000 /xm Lg-ioooMm > lO^^L©) were discovered to be a significant class of objects 
in the local Universe in the IRAS all sky survey (Soifer et al. 1987, Sanders et al. 1988). 
These systems have recently been suggested as a major component in the Universe at z> 2 
through deep submillimeter surveys ( c.f. Lilly et al. 1999, Blain et al. 1999, Barger et al. 
1998). 

Whether ULIRGs are powered by dust enshrouded AGN or starbursts has been debated 
since their discovery(e.g. Sanders 1999, Joseph 1999). This question is made all the more 
important by the realization that ULIRGs could produce a major fraction of the radiant 
energy in the Universe (Lilly et al. 1999, Heckmann 1998). 

A wide variety of observational techniques have been employed to probe the nature of 
the central luminosity source in ULIRGS, ranging from optical spectroscopy (Sanders et 
al. 1988) to cm radio continuum imaging (Condon et al. 1991, hereafter C91), to VLBI 
observations (Smith, Lonsdale and Lonsdale 1998) and hard X-ray observations (Rieke 
1988, Nakagawa et al. 1998, Isawawa 1999). The nuclear environments of these galaxies 
have been effectively imaged with high spatial resolution in the optical and near infrared 
(Surace et al. 1998, Surace and Sanders 1999, Scoville et al. 1999, hereafter S99). 

Observations in the thermal infrared have the potential to address the question of the 
origin of the luminosity by establishing the size of the regions that arc radiating the infrared 
luminosity. The IRAS observations showed that the mid and far infrared luminosity in these 
systems emerging predominantly between 50 and 200 /xm was due to thermal reradiation 
by dust. However, the spatial resolution of those observations, ~ 1' — 2', was inadequate 
to place meaningful size constraints on the luminosity sources. Because the mid-infrared 
wavelengths (8-25//m) carry a significant fraction of the total bolometric luminosity in 
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ULIRGs, ranging from ~10% to > 30% of the 8-lOOOyum luminosity, observations at these 
wavelengths with a resolution < 1" have the potential to constrain the sizes of the emitting 
regions and thereby address the nature of the underlying sources that heat the radiating 
dust. Observations in the 10/um atmospheric window, with a resolution ~ 1", have been 
reported for several ULIRGs by Miles et al. (1996), and Keto et al. (1992). In most cases 
they found the sources to be unresolved implying that the sources were either very compact 
star bursts or dust enshrouded AGN. 

With a diffraction limit of 0.24" FWHM at 10//m, the Keck Telescope provides a 
substantial improvement in spatial resolution over previous mid-infrared observations, 
probing the distribution of the thermal emission at the 100-300 pc scale in nearby ULIRGs. 
With the introduction of mid-infrared imaging on the Keck Telescope, we have started a 
program to address this problem. In the first paper from this program we reported imaging 
of the nucleus of the closest ULIRG - Arp 220 (Soifer et al. 1999). The observations were 
able to spatially resolve the emission at 24.5/um, and demonstrated that the central thermal 
source is optically thick at this wavelength. These observations placed severe constraints on 
the surface brightness in the infrared luminous nuclei. 

In this paper we report imaging observations from 8 - 24.5/i.m of a sample of the 
closest ULIRGs at spatial resolutions of 0.3-0.6". These data provide the highest spatial 
resolution yet achieved on the thermal emission from the closest ULIRGs and trace the 
spatial distribution of the emergent luminosity in these systems. We adopt Hq —75 km 
s-^Mpc-^ 
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2. The Sample 

The objects observed were taken from the IRAS Bright Galaxy Sample (Soifer et 
al. 1987) with the addition of IRAS 17208-0014 (Sanders et al. 1995), which is at a low 
Galactic latitude. The basic information for the seven objects in the sample is given in 
Table 1. The objects selected for observation were chosen on the basis of 1) meeting the 
luminosity definition of the ULIRGs, i.e. Ls-wooum > lO^^i^©, 2) having a 25//m flux 
density > 1 Jy based on IRAS measurements, 3) being at the lowest possible redshift for the 
best linear resolution, and 4) being available at the time of the observations. This sample 
contains the five closest ULIRGs in the IRAS Bright Galaxy sample of Soifer et al., and 
seven of the ten ULIRGs known with 5 > —30 deg and z < 0.06. The projected hnear scales 
of these objects range from 360 pc arcsec"^ to 1200 pc arcsec"^, so that the point spread 
function/ angular resolution of the Keck observations, 0.3" at 12//m and 0.6" at 24.5/im 
provide linear resolutions of 100 pc - 700 pc, depending on the object and the wavelength 
observed. Three of the objects in the sample have been classified as Seyfert nuclei, three as 
LINERS, and one as an HII region, based on optical spectroscopy. 

3. Observations and Data Reduction 

The observations were made primarily using the MIRLIN mid-infrared camera (Ressler 
et al. 1994) at the f/40 bent Cassegrain visitor port of the Keck II Telescope. The camera 
uses a 128x128 Si: As array with a plate scale of 0.138"/pixel for a total field of view of 
17"xl7". Table 2 reports the central wavelengths for the filters through which each object 
was observed. At each wavelength the observing procedure was the same. A secondary 
with a square wave chop of amplitude 6" in the north-south or east-west direction at 4 
Hz was employed for fast beam switching. The frames sampling each chop position were 
coadded separately in hardware, resulting in two images. After an interval of approximately 
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a minute, the telescope was nodded perpendicular to the chop direction (east-west or 
north-south) by 6" and a second pair of images was obtained in order to cancel residuals 
in the sky and to subtract telescope emission. This procedure was repeated a number of 
times at each wavelength. The data were reduced by differencing the two images obtained 
within the chop pairs at each nod location, and then coadding the resulting positive images, 
with the positions appropriately adjusted to a common location, to yield a positive image 
centered in a field approximately 6"x6". Because of the chopper and telescope nod spacings 
employed for the observations, the data are not capable of measuring flux outside a 6" 
diameter region. 

Observations with MIRLIN of all the ULIRGs except IRAS 05189-2524 were obtained 
in two nights in March 1998. The observations of the targets were interleaved with 
observations of nearby bright stars that were used to establish the point spread function 
(PSF) for the observations. The PSF stars were observed only three times per night in 
each of the bands at 12.5, 17.9 and 24.5 yum, but an attempt was made to accompany the 
critical observations with measurements of a nearby PSF star. At 24.5 fim the PSF was 
always equal to the diffraction limit of the telescope to be 0.62" FWHM. At the shorter 
wavelengths, the size of the PSF image was set by diffraction and the pixel sampling about 
half of the time; at other times atmospheric seeing affected the image size. 

IRAS 05189-2524 was observed in October 1998. During this run the PSF stars 
were observed more frequently, but for wavelengths shorter than A ~ 17.9 fxm significant 
variation in the measured size of the PSF star was seen. 

The variation in the measured size of the PSF calibration star is the largest uncertainty 
in the measurement of the source sizes; this issue is addressed when observations of a 
particular object are described. 

The MIRLIN observations were made under photometric conditions. The photometry 
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was calibrated based on observations of four bright stars, a Tau, a Boo, /? Peg, and a Get 
whose magnitudes, in turn, were based on IRAS photometry. The uncertainties in the 
photometry, based on the internal consistency of the observations, is estimated to be 5% at 
A <17.9 //m and 10% at 24.5 /im. The flux density corresponding to 0.0 mag (Vega-based) 
was taken to follow the prescription given in the Explanatory Supplement to the IRAS 
Catalogs and Atlases (Joint IRAS Science Team, 1989), and is given in Soifer et al. (1999). 

In addition to the MIRLIN observations, UGC 5101, Markarian 231 and 
IRAS 08572+3915 were observed using the Long Wavelength Spectrograph on the 
Keck I Telescope on the night of 01 April 1999. The night was not photometric and so 
only measurements of the sizes of the compact sources were made; all the observations were 
carried out at 12.5 iim. The chopper was set to an amplitude of 5" chopping north-south 
at a frequency ~ 5 Hz. Observations were made in a fashion similar to the MIRLIN 
observations, except that the telescope nodding was in the same direction as the chopping 
and the nodding amplitude was the same as the chopping amplitude. The data were 
reduced in a manner similar to the MIRLIN data. The pixel scale is 0.08" pixel"^. 

Observations of Markarian 273 were made at 3.45 /xm using the near-infrared camera 
at the Cassegrain focus of the 200-inch Hale Telescope at Palomar Observatory in May 
1999. The camera was read out in a 64 x 64 pixel subarray because of the high read rate 
required to avoid saturation. The scale is 0.125" pixel"^. The chopper throw was 15" at a 
frequency of ~1 Hz. Observations of stars given by Elias et al. (1982) were used to estabhsh 
the photometric scale. 
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4. Results 

The mid-infrared observations of the ULIRGs reported here show a variety of 
characteristics, generally consistent with their properties seen at high spatial resolution at 
radio wavelengths (c.f. C91) and in the near infrared using NICMOS on HST (Scoville et 
al. 1998, S99). The mid-infrared data on all the ULIRGs are presented in Figures 1, 2 and 
3. In Figure 1 we present the 12.5 //m images of the ULIRGs and, where the radio and near 
infrared data arc available, compare these maps to the equivalent maps at 8.4 GHz from 
the data of C91 and to the NICMOS 2.2 /im data from S99 and Scoville et al (1998). In the 
case of three galaxies, IRAS 05189-2524, IRAS 08572+3915 and Markarian 231, where there 
is a single compact nucleus that emits a substantial fraction of the 12.5/xm flux density, 
we have presented the corresponding PSF star at 12.5//m for comparison. In Figure 2 we 
compare, for the five single nucleus sources, the integrated flux at 12.5 /xm as a function 
of beam radius with the same quantity for a corresponding point source to establish the 
spatial distribution of the mid-infrared radiation. Figure 3 displays the spectral energy 
distributions (SEDs) obtained from the mid-infrared data and compares these data with 
the corresponding IRAS data and data at other wavelengths. The flux densities measured 
for the ULIRGs in a 4" diameter beam are presented in Table 3. 

The mid-infrared observations did not estabhsh accurate absolute positions for the 
sources. In some cases, described below, we use arguments regarding the similarities of the 
structures seen at vastly different wavelengths to argue for the physical coincidence of the 
sources seen at different wavelengths. 

In the following sections we give results for the individual objects. 
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4.1. IRAS 05189-2524 

MIRLIN observations of IRAS 05189-2524 were made at 12.5 and 24.5 pim. The 
shorter wavelength data have the higher signal-to-noise ratio, and so are presented in Figure 
1. Images at 12.5/^m (MIRLIN), 2.2/im (NICMOS, S99) and 8.4GHz (VLA, C91) are shown 
in Figure la. The dominant central source is clear in all three maps. At 8.4GHz, the single 
source is concentrated with a measured size of 0.2"x0.17"(C91). At 2.2 //m, the unresolved 
source is < 0.2" and contains ~75% of the total flux measured in a 5" diameter beam, and 
70% of the flux contained in an 11" diameter beam (S99). 

The source appears virtually unresolved in the MIRLIN image. As can be seen from 
Figure 2, there is excellent agreement between the curve of growth for IRAS 05189-2524 
and the PSF star observation. The apparent size of the PSF star varied during the time of 
the observations, so the hmit on the size of IRAS 05189-2524 is not as stringent as would 
otherwise be achieved. At 12.5/xm the size is <0.2" (Table 4). The source size at 24.5 //m is 
an upper limit as well, and is larger and thus consistent with that established at 12.5 fj,m. 

As seen in Figure 3 and Table 3, the flux densities measured by MIRLIN at 12.5 /^m 
and 24.5 yum agree well with the IRAS observations of this source; virtually all of the 
luminosity of this source at 12.5 fj,m and 24.5 fj,m is contained within the unresolved source 
measured here. 

The MIRLIN observations of IRAS 05189-2524 lead to a limit on the source diameter 
of <170 pc in the mid-infrared (Table 5). The combination of flux density and upper limit 
on the source sizes(Table 4), imply the brightness temperatures T^ >123 K and >85 K at 
12.5 /im and 24.5 fj,m respectively. The flux ratio between 12.5 /im and 24.5 fj,m corresponds 
to a color temperature Tg =169 K, while the color temperature derived from the 25/um and 
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60//m IRAS measurements is 83 K^. 

A central source with the bolometric luminosity of IRAS 05189—2524 would heat 
silicate or graphite grains to the apparent color temperature between 12 and 24.5//m, Tc ~ 
170 K, at a radius of 30-60 pc using the grain parameters of Draine and Lee (1984). This is 
much smaller than the observational limit. Alternatively, the size of an optically thin source 
which produces the 25—60 fxm flux ratio with a color temperature of 83 K is significantly 
larger, 170-250 pc in radius, again using the Draine and Lee grain parameters. 



4.2. IRAS 08572+3915 

MIRLIN observations of IRAS 08572-^3915 were obtained at five wavelengths (Table 
2). In addition, LWS observations of IRAS 08572-1-3915 were obtained at 12.5 /im. Optical 
and near infrared imaging (Sanders et al. 1988, Carico et al. 1990, S99) show two galaxies 
separated by about 5" in a northwest-southeast orientation. The southeast source is not 
detected in the radio (C91) or in the present data. 

The 12.5/xm LWS image in Figure lb shows a bright central source (the northwest 
source in the near infrared) that is unresolved. The size limit at 17.9 fim is consistent with 
that found at 12.5 //m. Figure lb and Figure 2 show that the central source is dominant 
in the mid-infrared and radio, containing ~80% of the total flux measured in a 4" diameter 
beam. At 2.2 ^m, the central unresolved source contains ~60% of the total flux in an 11.4" 
diameter beam (S99). 

^The brightness temperature, T;,, is defined hy fi, = fl x B,^{Ti,) where B^^iTj^) is the 
blackbody emission at frequency v for temperature T^ and is the solid angle subtended 
by the source. The color temperature of a source, Tc, measured between frequencies vi and 
V2 is found from 4^ = ^i^^4Sv. 



- 12 - 



The most stringent observations constraining the size and spatial distribution of the 
source are the LWS observations at 12.5/xm. Five observations of the PSF star (BS3275) 
were interleaved with eight observations of the object. The average FWHM of the object 
(0.39 ± 0.04") was marginally greater than the average FWHM of the PSF (0.34 ± 0.01"). 
Formally this leads to a size of the 12.5//m source of 0.19", but this is only a lu result. A 3a 
limit on the size of the compact source is ~0.3". As illustrated in Figure 2, there appears to 
be a measureable flux, approximately 20% of the total at 12.5 /im, that extends beyond the 
central core, and is contained within 4". The figure shows clearly that the fiux measured in 
IRAS 08572+3915 continues to increase beyond the beam where the flux in the PSF star 
has stopped increasing. 

Figure 3 and Table 3 compare the Keck photometry in a 4" diameter beam for 
IRAS 08572+3915 with the IRAS measurements. The MIRLIN photometry accounts for 
70% of the IRAS flux density at 12.5y[/m. The MIRLIN observations were taken with 
comparatively narrow fllters, i.e., AA/A <~ 0.1, from 7.9yum to 17.9//m, while the IRAS 
12/im filter bandpass extends from 8-15^m. Given the uncertainties in the corrections for 
the different bandpasses, the MIRLIN data are consistent with all the fiux detected by 
IRAS being confined to a 4" diameter beam surrounding the detected source although it is 
possible that some of the fiux could be accounted for by extended emission associated with 
the northwest source or from the southeastern source. 

Neither the MIRLIN nor the LWS observations detected the galaxy approximately 5" 
southeast of the object centered in Figure lb. The MIRLIN data limit the contribution of 
the southeastern source to <30% of the total 12.5yum fiux; it is hkely that the southeast 
galaxy does not contribute significantly to the total infrared luminosity in this system. This 
is consistent with the optical and near infrared morphology and colors of this galaxy, which 
are indicative of quiescent galaxy nucleus (Surace et al. 1998, Surace and Sanders 1999). 
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4.3. UGC 5101 = IRAS 09321+6134 

MIRLIN observations of UGC 5101 were obtained at 12.5 and 17.9 //m. Contour plots 
of this source in Figure Ic show that at all wavelengths the structure is basically the same: 
a bright central core centered on an apparently elliptical structure having a major axis at 
a position angle of ~80°. While the position angle of the elliptical structures agree well at 
2.2yum and 8.4 Ghz, there is a suggestion that the position angle of the ellipse at 12.5 //m is 
~15° smaller. At 8.4 Ghz the central source contains 80% of the radio fiux(C91), while at 
2.2//m, the unresolved peak contains 35% of the flux in a 5" diameter beam (S99). 

As shown in Figure 2, at 12.5 fim approximately half of the total flux measured in a 
4" diameter beam is contained in the central compact source. The inner central compact 
source is indistinguishable from the PSF (Figure 2). From the dispersion in the sizes 
determined for the PSF star (BS3775) the compact source has a size (3cr) less that 0.22". 
The size measured at 17.9//m is consistent with that measured at 12.5 //m. 

The data of Figure 3 and Table 3 show that at 12.5 nm the flux measured by MIRLIN 
is ~60% of that measured by IRAS. This suggests that there is low surface brightness 
emission extended on a larger scale than 4". The 12.5fj,m flux density contained in the 
compact source is ~30% of the total flux density at 12//m measured by IRAS. 

4.4. Markarian 231 = UGC 8058 = IRAS 12541+5708 

Markarian 231 has long been known both as an infrared luminous system (Rieke and 
Low, 1972) and a Seyfert 1 system (Boksenberg et al. 1977). The IRAS survey found that 
it is the most luminous object within 300 Mpc, with a bolometric luminosity of 3 x lO^^L©. 

MIRLIN images of Markarian 231 were obtained at seven wavelengths from 7.9 to 
17.9//m. In addition, LWS imaging was obtained at 12.5/xm. Figure Id shows contour plots 
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of Markarian 231 at 12.5 //m from the LWS data and data at 8.4 Ghz(C91). The emission 
at all the mid-infrared wavelengths is dominated by the very bright central source, and is 
effectively unresolved. The radio continuum is entirely from the bright nucleus. 

LWS imaging was used to establish the size of the central source at 12.5/xm. The LWS 
observations consisted of eleven measurements of the size of Markarian 231, interleaved with 
nine measurements of a near-by PSF star(BS 4905) over a period of two hours. The average 
FWHM of the PSF star and Markarian 231 were indistinguishable and were shghtly greater 
than the diffraction limit of the telescope and the rms (population) scatter in the FWHM 
was 0.04". The central source thus is unresolved in these observations. We have chosen to 
set a 3(7 limit on its FWHM of 0.13" on the basis that a difference three times the average 
standard deviation of the FWHM would have been easily detected. 

The SED of Markarian 231 is plotted in Figure 3, and shows clearly that all the flux 
measured by IRAS is detected in the 4" diameter beam. The SED shows a dip at 10//m 
apparently due to foreground silicate absorption of the central source. 

4.5. Markarian 273 = UGC 8696 = IRAS 13428+5608 

The MIRLIN observations of Markarian 273 consist of imaging at seven wavelengths 
from 7.9//m to 17.9//m. Figure le shows a comparision of the MIRLIN 12.5//m data with 
the 8.4 GHz map (C91) and the NICMOS 2.2/im image (S99), while a montage of images of 
Markarian 273 is shown in Figure 4. At 12.5/im the source structure is basically double, with 
the northern source brighter than the southwestern source at most wavelengths. A similar 
morphology is seen at 2.2/xm (S99), but not in the radio (C91). In comparing the 12.5//m, 
NICMOS and 8.4 Ghz data we have assumed the spatial coincidence between the bright 
northern source in each image. The overall ageement between the 12.5/im and NICMOS 
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images is excellent, while the discrepancy between the positions of the southwestern source 
in the infared and radio is real. The position angle between the northern and southwestern 
sources is the same (to within the uncertainties) in the 12.5/im, 2.2//m and 8.4 GHz 
images, while the source separations are 1.0", 1.0" and 1.2" at these same wavelengths. The 
apparent separation between the infrared and radio sources, 0.2", corresponds to 160pc at 
the distance of Markarian 273. Note also that the flux ratio between the northern and 
southwestern source is considerably different between the infrared and radio. 

As shown in Figure 4, the structure seen in Markarian 273 at all infrared wavelengths 
from SAfim to 17.9/Ltm is basically the same. There is no evidence for a mid-infrared 
counterpart to the southeastern source seen at 8.4 Ghz by C91 (Figure le) and at 5 Ghz 
by Knapen et al. (1997). The bright source present to the southeast in the radio map 
is suggested to be a background radio object by Knapen et al., although the statistical 
hkelihood of a background source faUing within 1" of the location of the radio core of 
Markarian 273 is < 10~^. There is a faint compact object found in the NICMOS images 
(S99) at the location of southeast radio source, although there is no evidence for any 
mid-infrared emission at this location. The upper limit on the 12.5/xm flux density at the 
location of the southeastern radio source of 15 mJy is 1/10 that of the northern source. 
This hmit, coupled with the 8.4GHz flux density of this source of 3.2 mJy (C91), leads to 
a limit on the ratio of flux densities ^"(840^) < 4.7, consistent with the observed ratio 
f^is^AOHz) =4.5 for the northern source. By comparison, if the radio ( 0.2mJy at 8.4 GHz, 
C91, Knapen et al.) and infrared (58 mJy at 12.5//m) flux densities of the southwestern 
source are assumed to be associated with the same physical source, the ratio y''(84Gjj^) — 
290, a factor of 65 times greater than that seen in the northern source, making this nucleus 
far more infrared bright relative to the radio than the northern source. If the southwestern 
radio and infrared sources are not related, the southwestern infrared source has an even 
larger ^§5^ ^^tio. 
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As can be seen in Figures le and 4, the relative brightnesses of the nuclei change 
substantially as a function of wavelength. In order to get an estimate of the relative flux 
emitted by the two nuclei, the flux in two identically sized square beams was measured at 
the same locations on the images at all wavelengths presented; the beams, 1.6" on a side, 
are shown in Figure 4. The sky was taken northwest and southeast of each of the squares. 
The efficiency of the square beams was estimated by placing the PSF of the appropriate 
wavelength at the locations of the nuclei within the square beam. The photometry of the 
NICMOS image was arbitrarily normalized so that the sum of the fluxes from the two 
squares equals the sum of the fluxes give by Knapen et al. (1997) and the flux density 
in a 4" diameter beam equals that of Knapen et al. The results of these photometric 
measurments are shown in Figure 5. 

At 3.4//m and 10.3//m the southwestern source is slightly (10%) brighter than the 
northern source. At the other mid-infrared wavelengths the northern source is between a 
factor of 1.5 and 3 times brighter than the southwestern source. The different SED's of 
the components implied by their changing relative brightnesses is shown directly in Figure 
5, where photometry of the two nuclei is shown. The northern source clearly shows a 
signflcantly deeper absorption at 10//m than the southwestern source. This is consistent 
with the distribution of CO emission which strongly peaks on the northern source (Downes 
and Solomon, 1998), but is inconsistent with the near infrared colors found by S99, where 
the northern source has a smaller J-K color than does the southwest source. 

The SED of Figure 5, particularly the relatively strong emission at 3.5 and 11.7 /xm, 
suggests aromatic hydrocarbon emission (PAH) emission associated with the southwest 
nucleus. Such emission is seen in the spectra of many ULIRGs (Genzel et al. 1998). The 
emission from the southwest source could account for all of the 11.3 //m emission seen in 
the spectrum of the combined nuclei by Dudley(1999).The northern source does not show a 
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similar peak at 11.7yum, but clearly shows the effects of substantial overlaying cold silicate 
absorption. 

The MIRLIN photometry shown in Figure 3 was obtained using a 4" diameter beam 
centered between the two nuclei. The comparison of the MIRLIN and IRAS data (Figure 
3) demonstrates that the mid-infrared emission in this source emerges virtually entirely 
from the central region. Because of the proximity of the two nuclei and the fairly low 
signal-to-noise ratio of the nuclei in some of the images, it is difficult to cleanly determine 
the sizes of the two sources. The data suggest that both nuclei have unresolved mid-infrared 
cores with sizes <0.4". In the northern nucleus approximately 2/3 of the 12.5//m flux is 
unresolved, while in the southwestern nucleus, about half the flux is unresolved. 

4.6. Arp 220 = UGC 9913 = IRAS 15328-|-2340 

The data for this source were discussed in detail in Soifer et al. (1999) and are included 
in Figures If and 3 for completeness. 

4.7. IRAS 17208-0014 

MIRLIN imaging was obtained at 7.9 and 12.5/im. The 12.5/im image is shown in 
Figure Ig. The 2.2//m NICMOS image of IRAS 17208-0014 (Figure Ig, S99) shows it to be 
an extended source of roughly 2" diameter with a number of luminous starburst centers to 
the northeast of the nucleus. The mid-infrared image is apparently also similar, but the 
limitation due to the chop-nod spacing and the relatively low infrared surface brightness 
results in an image with low signal-to-noise ratio. The size of the 12.5 fj,m source of ~2" 
or 1.6 kpc FWHM (Figure 2) is quite similar as well to that found in CO by Downs and 
Solomon (1998) and Jogee et al.(2000). 
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The best illustration of the spatial distribution of the mid-infrared emission is shown 
in Figure 2. The source flux smoothly increases up to 5" diameter, with a FWHM of 2"; 
the apparent leveling off of the enclosed flux at a diamteter of 5" is caused by the artificial 
imposition of the sky level between diameters of 5 and 6". 

Evidence that the source is not extended much beyond a diameter of 5" is provided by 
Figure 3 which shows that, within the uncertainties, the enclosed flux within a 4" diameter 
beam is the same as within the IRAS beam. The lack of a dominant compact component 
in the mid-infrared make this object very different from the other ULIRGs in this sample. 

The radio emission from IRAS 17208-0014 is much more compact than the infrared 
emission. At least 90% of the 1.4GHz flux comes from a region smaller than 2"x2" (FWHM, 
Condon et al. 1996). Martin et al. (1989) measure a radio source with a size of 0.26"x0.32" 
at 1.6 GHz that accounts for ~90% of the 1.4 GHz flux density, suggesting a compact radio 
core. This core is clearly not dominant in the mid-infrared imaging. The near infrared 
NICMOS imaging of this galaxy (S99) suggests that the nucleus could be obscured at 
2.2/xm, but the lack of detection of a strong compact core at 12.5//m is quite peculiar, since 
the radio source is so compact. 

5. Discussion 

5.1. Compact Sources in ULIRGs 

The data presented here, with angular resolution of <0.3", represent the highest spatial 
resolution achieved to date in probing the central regions of ULIRGs in the thermal infrared. 
Table 4 summarizes the key observations, the sizes of the compact sources found in the 
ULIRGs at 12.5 //m, 17.9 /xm and 24.5 //m, as well as the fraction of the total flux at these 
wavelengths emitted from these cores. All of the ULIRGs except IRAS 17208—0014 have 
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a nuclear source in the mid-infrared of size <200 pc that produces a substantial fraction, 
if not all, of the galaxy's mid-infrared luminosity. In three sources, IRAS 05189—2524, 
IRAS 08572-1-3915 and Markarian 231, a single unresolved nucleus produces between 80% 
and 100% of the total mid-infrared luminosity. In Markarian 273 and Arp 220, two nuclei 
produce nearly all the mid-infrared emission in the galaxies. In the case of Markarian 273 
the nuclei are substantially unresolved, while in Arp 220 at least one nucleus is resolved. 
Because of its low redshift the linear size of the resolved nucleus in Arp 220 is smaller 
than the limit on all the unresolved nuclei except Markarian 231(see Table 5). In only one 
compact system, UGC 5101, does the compact nucleus produce significantly less than half 
the total mid-infrared emission from the system. 

As noted previously, IRAS 17208—0014 presents an interesting puzzle. There is no 
mid-infrared evidence for a compact nuclear source, but the radio data (C91, Martin et al. 
1989) strongly suggest a power source with a size <0.3". In the following discussion, we 
will adopt the morphology implied by the 12.5//m image, but we note that extremely high 
12.5yum extinction could be leading to an erroneous interpretation of the nuclear structure 
in IRAS 17208-0014. 

Table 5 shows that the compact nuclear sources have linear sizes typically <200 pc in 
diameter. In Arp 220 the diameter of the western source is 140 pc at 24.5//m (Soifer et al. 
1999), while in Markarian 231, the diameter of the 12.5/xm source is <100 pc. 

5.2. AGNs as the Luminosity Source in ULIRGs 

One of the goals of the observations reported here is to distinguish the ultimate source 

of the infrared luminosity in these systems by means of the spatial extent of the infrared 
emission. If the dust is heated by a bright compact object such as an AGN, the dust 
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emission will be compact, and the size of the infrared source will increase with increasing 
wavelength roughly as r ~ A^-^, since the dust temperature will decrease with radius as 
T ~ Qj^^ lY^Q wavelength of maximum emission X^ax ~ T~^. If, on the other hand, 
the heating is due to many hot stars in a starburst region, then the dust emission will be 
distributed on the same scale as the young, luminous stars. Since the dust heating is due to 
the proximity of the dust to the luminous young stars there will be little or no temperature 
gradient seen. Spatially extended 10 /xm emission therefore implies starburst activity, while 
spatially compact 10 //m emission could be due to either a compact starburst or an AGN. 

To quantify the compact mid-infrared emission we utilize the brightness temperatures 
derived from our observations and the color temperatures determined from the IRAS 25 
and 60/xm observations of these objects. We assume that Tc measures the temperature of 
the grains heated by the luminosity source; should be <Tc if the optical depth t^, < 1. 
Only in the limit of large optical depth should the brightness temperature approach the 
color temperature. 

Based on the MIRLIN observations of Arp 220, Soifer et al. (1999) showed that the 
dominant luminosity source, Arp 220W, is optically thick with a brightness temperature at 
25 /um equal to the color temperature between 25 and 60 /xm. The observations reported 
here permit analogous constraints to be placed on some of the nuclear sources in other 
ULIRGs in this sample. In Markarian 231, IRAS 05189-2524 and IRAS 08572+3915, the 
limits on the sizes of the central infrared sources are sufficiently small that it is worthwhile 
considering how these sizes constrain the far infrared sources. 

The size limit on Markarian 231 of 0.13" at 12.5^m places the most severe limits on 
models of the infrared emitting region. In this case the observed flux density and size 
hmit lead to a lower hmit on the brightness temperature of T;, > 141K at 12.5//m. If the 
infrared emitting source has the same size at 17.9 fim, 25 fj,m and 60 /im, the corresponding 
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brightness temperatures are 124 K, 128 K and 241 K. The brightness temperatures at 
12.5/xm, 17.9//m and 25 /im are consistent with each other and with a nearly optically thick 
source over this wavelength range, but the 60 /im brightness temperature is much too large 
to be consistent with dust models with reasonable grain properties. This implies that the 
60/xm source size is much larger than the 12.5//m size. 

The color temperature measured between 25 and 60 //m, determined from the IRAS 
measurements, is 85 K. The size of a blackbody of 85 K required to produce the observed 
60/im flux density is 0.4", significantly larger than the limit on the source size at 12.5/um. 
The source size determined if Tc = T;, represents a lower limit on the size of the 60/xm 
source. For example, if the size of the 60//m source were equal to the ~1" size found in 
CO by Downs and Solomon (1998), the 60yum brightness temperature would be 53K, much 
less than that inferred by assuming all the 25 and 60 /um emission is emerging from the 
same region. The size limits are based on the assumption of large optical depth (r >>1); 
smaller optical depth would imply an even larger source. Since significant emission is seen 
at A < 60/xm, either reo^m <1 or we are viewing a tilted optically thick disk, as suggested 
for Arp 220 (Scoville et al. 1998, Downs and Solomon, 1998, Sakamoto et al. 1999, Soifer 
et al. 1999). In either case, the size constraints show that in Markarian 231 the infrared 
source increases in size with increasing wavelength. This is consistent with the picture of 
an AGN with luminosity of 3x10^^ Lq heating the surrounding dust and producing a radial 
temperature gradient. 

In this sample there are two additional ULIRGs where significantly more than 50% of 
the mid-infrared luminosity emerges from a compact central source, IRAS 05189-2524 and 
IRAS 08572+3915. In IRAS 05189-2524, if the size of the 60//m emitting source is the same 
as the limit on the 12.5/im source measured here, Tfe(60//m) = 101 K, which is significantly 
greater than Tc =85 K from the IRAS 25 and 60/Ltm data. For the same assumptions, in 
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IRAS 08572+3915 T6(60/xm)=77 K is slightly less than Tc =83 K determined from the 
IRAS data at 25 and 60/im. For both sources, using dust properties described by Draine 
and Lee (1984) the sizes calculated for silicate and graphite dust grains heated by central 
sources of the observed bolometric luminosity, and having dust temperatures equal to the 
observed color temperatures are larger by a factor of 1.5 to 3 than the observed sizes of 
these sources at 12.5//m. This suggests that the sizes of these infrared sources also increase 
with increasing mid and far-infrared wavelength. 

High apparent surface brightnesses and a radial temperature gradient are natural 
consequences of centrally heated dust emission surrounding a high luminosity source, i.e., 
an AGN. In Markarian 231 there is direct evidence that the size of the infrared source 
increases with increasing wavelength in the mid and far infrared. In IRAS 05189-2524 and 
IRAS 08572-1-3915, there is strong evidence for such radial temperature gradients based on 
the compact sources observed in the mid-infrared. The observed luminosities are adequate 
to quantitatively account for the inferred variation of dust temperature with radius if the 
dust is heated by a central source in these objects. 

5.3. Starbursts as the Luminosity Source in ULIRGs and Energy Production 

Rates 

The compact nature of the sources in the objects observed leads to very high 
luminosities produced in small volumes. By comparing these energy production rates with 
those found in local examples of starburst environments, we can determine whether known 
starburst environments can generate the apparent energy production in the ULIRGs. 

In the core of Orion the luminosity of 2 x lO^L© emerges from a 0.3 pc size region 
(Werner et al. 1976), while in the nuclear starburst in M82, the nearest starburst galaxy, the 
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luminosity of 3 x lO^^L© is emitted over a projected area of 450x75pc (30" x 5") (Rieke et 
al. 1980). These observations lead to apparent surface brightnesses of emergent luminosity 
of 1 — 3 X 10^^ L0/kpc^. If, as seems likely, M82 is being viewed edge-on, then the mean 
surface brightness of this system, if viewed face on, is closer to 2x10^^ Lo/kpc^. 

A survey of the maximum surface brightness in nearby starburst galaxies in the UV, 
infrared and radio by Meurer et al. (1997) leads to peak global averaged surface brightnesses 
in starburst galaxies of ~ 2 x 10^^ Le/kpc^. Meurer et al. find that the peak surface 
brightnesses of clusters in starburst galaxies reaches ~ 5 x lO^^LQ/kpc^ over regions ~10 
pc in diameter. Meurer et al. suggest that the global peak surface brightnesses is due to 
a mechanism whereby the pressure from stellar winds and supernovae is able to regulate 
the star formation rates in galaxy disks, while the formation of massive stars in giant HIT 
regions is regulated by the "Kennicutt Law" , i.e. star formation occurs when the surface 
gas density exceeds a critical value that depends on the rotational velocity and shear 
(Kennicutt, 1989). 

The surface brightnesses observed in the ULIRGs can be estimated if we assume 
that the size (or limit) of the region measured in the mid-infrared is representative of the 
infrared emitting region as a whole. With this assumption, wc have presented the surface 
brightnesses for the observed ULIRGs in Table 5. We include two values of the surface 
brightness, a "maximum" which assumes that the fraction of the total infrared luminosity 
emerging from the compact core is the same as the fraction of the mid-infrared flux density 
emerging from the core, and a "minimum" that assumes that only the observed mid-infrared 
luminosity emerges from the compact core. In the case of Markarian 231, the former 
number is unphysically large (section 5.2). In the other systems, the two values represent 
likely upper and very conservative lower limits on the apparent surface brightness of the 
infrared emission. 
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Only in the case of IRAS 17208—0014, where the near and mid-infrared observations 
suggest the luminosity is due to star formation, is the measured infrared surface brightness 
comparable to the peak global value 2 x 10^^ LQ/kpc^ found by Meurer et al.(1997). Other 
than this case, the inferred surface brightnesses are much greater than found globally in disk 
starburst systems. The highest apparent surface brightnesses in the ULIRGs are found in 
the western nucleus of Arp 220, which is most likely powered by a nuclear starburst (Smith 
et al. 1998b) and the cores of Markarian 231 and IRAS 05189—2524, which are probably 
powered by AGN. In these cores the apparent surface brightness approaches or surpasses 
the surface brightness of clusters in starburst galaxies, while in the other systems (excluding 
IRAS 17208—0014) the surface brightness most probably exceeds 10^^ L0/kpc^. In all these 
cases the surface brightness over regions of hundreds of parsecs across are comparable to 
the brightness found in ~10 pc sized clusters within starburst systems by Meurer et al. 

For the ULIRGs that are powered by star formation, the large surface brightnesses 
and sizes suggest that the global starburst intensity regulating mechanism of Meurer et al. 
(1997) does not apply. Gravitational instability allows for higher gas densities and hence 
higher star-formation rates in the solid-body (no shear) rotation nuclear regions. Since the 
nuclei of most ULIRGs are mergers, far from equilibrium, the gas densities may be very 
high, leading to high star formation rates and correspondingly high surface brightnesses for 
short times. ULIRGs are clearly extraordinary events in the evolution of galaxies. Several 
studies (e.g., Sanders et al. 1988, Scoville, Yun and Bryant 1997, Downes and Solomon 
1998) suggest that the nuclear disks are unstable in some ULIRGs. The large gas masses 
in these environments might then lead to nuclear star formation rates much greater than 
those found in less extreme systems. 
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5.4. A Panchromatic View of ULIRGs 

Table 6 summarizes the characterization of the ULIRGs applying many different 
diagnositics of the underlying energy sources. Here we have summarized the structure in the 
mid-ir from the data presented here, in the near infrared from S99 and Surace and Sanders 
(1999), and in the radio from the VLBI data of Smith et al. (1998). The optical spectral 
classifications are from Sanders et al. (1988) and Veilleux et al (1995). The near-infrared 
spectral classifications are from Veilleux et al. (1999) and Murphy (private communication). 
The ISO spectral classifications are based on the strength of the PAH features and are taken 
from Genzel et al. (1998) and Lutz, Veilleux and Genzel (1999). The X-ray classifications 
are from ASCA results and are taken from Iwasawa (1999) and Nakagawa et al. (1999). 

In none of the objects in this sample, excepting perhaps Markarian 273 and IRAS17208- 
0014, is there a completely consistent set of "classifications" based on the wide range of 
wavelengths over which these objects are observed. In the object that is the clearest case for 
a central AGN luminosity source, Markarian 231, the hard X-ray emission is not consistent 
with the simplest view of such power sources (Nakagawa et al. 1999). In the predominantly 
Seyfert nucleus IRAS 05189-2524 there is not a dominant compact radio core. In Markarian 
273 the two nuclei clearly seen in high resolution imaging complicates the classification, 
since none of the spectral classifications resolved the separate nuclei. In UGC 5101, another 
source where the optical, infrared and X-ray indications point towards a starburst power 
source, a compact radio core characteristic of a central AGN exists. Compact mid-infrared 
sources occur in six of seven of these sources, where other evidence points strongly to both 
AGN and starbursts as the dominant underlying power source. Table 6 shows that most 
ULIRGs do not present a picture over all wavelengths consistent exclusively with either 
an AGN or starburst power source. In most objects observed here, it is likely that both 
contribute to some extent to their total bolometric luminosity, though Markarian 231 and 
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IRAS 17208-0014 appear to represent the extremes of dominant AGN and starburst sources. 

6. Summary and Conclusions 

Diffraction limited mid-infrared observations of a sample of seven ULIRGs obtained on 

the Keck 10m telescopes have shown: 

1 - Extremely compact structures, corresponding to emission on spatial scales of 

< 100— 200 pc, are seen in six of the seven systems observed. The mid-infrared emission in 
these galaxies generally emerges from a single region of diameter < 1 kpc, and in most cases 
from regions <200 pc in diameter. 

2 - The upper limit on the diameter of the 12.5 //m source in Markarian 231, 0.13", is 
too small to be consistent with the observed 60 //m flux density emerging from a region 
this size, and demonstrates that in this AGN the source size increases with increasing 
wavelength (decreasing dust temperature). 

3 - In IRAS 05189-2524 and IRAS 08572+3915 there is strong circumstantial evidence 
that the size of the infrared source increases with wavelength between 12.5 and 60 //m, 
suggesting that these too are predominantly AGN powered ULIRGs. 

4 - In the other objects observed, the angular size limits for the compact mid-infrared 
sources are not yet adequate to distinguish between AGN and starburst luminosity sources, 
although the small sizes are consistent with a central luminosity source heating the 
surrounding dust as in AGN models in all the sources except IRAS 17208-0014. 

5 - If star formation powers the luminosity in the very compact infrared sources seen in 

these systems, the star formation rates averaged over a few hundred parsccs is comparable 
to that seen in the brightest ~10 pc sized clusters in nearby starbursts, and exceeds by 
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factors of up to 100, the global star formation rates in nearby starburst galaxies. 
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Fig. 1.- 

Figure la. Images of IRAS 05189-2524 obtained at 2.2/xm (S99, upper left), 8.4 GHz 
(C91, lower left), and 12.5yum (this work, upper right). The lower right panel is the PSF 
star image obtained at 12.5/xm contemporaneously with those of IRAS 05189-2524 and 
used to determine the 12.5//m size of IRAS 05189-2524. The peak of each image is at the 
location 0,0 in each panel. The contours in each panel are set so that the peak contour is 
90% of the maximum surface brightness. The contours are spaced down by a factor of 1.34 
so that the third contour down from the peak is at 50% of the peak surface brightness. 
The ninth contour down is 8.5% of the peak. Below this level the contours are spaced by 
a factor of two. The faintest contour is set at about 3 times the noise in the image. The 
mid-infrared data are from the MIRLIN observations. The hatched circle show the FWHM 
of the PSF for the 2.2/im observations. 

Figure lb. Images of IRAS 08572-1-3915. The format, source locations and contour 
levels are the same as in Figure la. The mid-infrared data are from the LWS observations. 
The 2.2//m and 8.4 GHz data are from S99 and C91 respectively. The hatched circles show 
the FWHM of the PSF for the appropriate observations. 

Figure Ic. Images of UGC 5101. The format is 2.2//m image(S99) - left, 12.5yumimage - 
center, 8.4 GHz image(C91) - right. The source locations and contours are set as in Figure 
la. The mid-infrared data are from the MIRLIN observations. The 2.2yum and 8.4 GHz 
data are from S99 and C91 respectively. The hatched circles show the FWHM of the PSF 
for the observations. 

Figure Id. Images of Markarian 231. The format is 12.5/im image - left, 12.5/xm PSF 
star - center, 8.4 GHz image - right. The source locations and contours are set as in Figure 
la. The mid-infrared data are from the LWS observations. The 8.4 GHz data are from C91. 
The hatched circle in the 8.4 GHz panel shows the FWHM of the PSF for this observation. 
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Figure le. Images of Markarian 273. The format, source locations and contour levels 
are the same as in Figure Ic. The mid-infrared data are from the MIRLIN observations. 
The 2.2//m and 8.4 GHz data are from S99 and C91 respectively. The + symbols are at 
the same coordinates in each panel, and are centered on the peaks of the 2.2//m image. 
There is excellent agreement between the locations of the 2.2//m and 12.5/im peaks. The 
discrepancy between the positions of the infrared peaks and the radio peak is real (see 
text). The hatched circles show the FWHM of the PSF for the observations. 

Figure If. Images of Arp 220. The format, source locations and contour levels are the 
same as in Figure Ic. The mid-infrared data are from the MIRLIN observations and are 
from Soifer et al. (1999). The 2.2//m and 8.4 GHz data are from S99 and C91 respectively. 
The hatched circle in the 12.5/xm panel shows the FWHM of the PSF for the observation. 

Figure Ig. Images of IRAS 17208-0014. The left image is 2.2/xm, the right image 
is 12.5/xm. The source locations and contour levels are the same as in Figure la. The 
mid-infrared data are from the MIRLIN observations. The 2.2//m data are from S99. The 
hatched circles show the FWHM of the PSF for the observations. 

Fig. 2.- 

The curves of growth for the objects (solid) and PSF's (dashed) for the five ULIRGS 
which are single nuclei. In all five the PSF's have been normalized so the enclosed flux of 
the PSF within a radius of 2.5" equals that of the object. For the objects, the sky level was 
set to be the value in an annulus between 2.5 and 3". In the plots for IRAS 08572-1-3915 
and UGC 5101 two normalizations are used for the PSF; in the lower curve the PSF is set 
to match the central profile of the galaxy, while the upper curve matches the enclosed flux 
from the PSF and the galaxy at the outer radius of the measurement beams, as for the 
other galaxies. The data of 1RAS0857+39 and Markarian 231 were obtained with LWS as 
outlined in the text; the other three data sets are from MIRLIN observations. 
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Fig. 3.- 

The spectral energy distributions of the seven ULIRGS studied are shown as a function 
of the rest wavelength from ~1 //m to 1 mm. Filled circles represent the MIRLIN data, 
while open circles have been taken from Carico et al. (1988) (near infrared data, 5" diameter 
beam) and Soifer et al. (1989) (IRAS data ~1' diameter beam). The MIRLIN photometry 
has been calculated for a 4" diameter beam centered on the object or between the two 
nuclei if multiple nuclei are present. As discussed in the text, the photometry is based on 
IRAS standards. 

Fig. 4.- 

Montage of images of Markarian 273. The images at wavelengths from 7.9//m through 
17.9/xm obtained with MIRLIN. The image at 3.4^m was obtained on the 200-inch Hale 
Telescope using the Cassegrain infrared camera. The images are centered so that the 
northern peak is at the location 0,0 in each panel. The locations of the boxes in the 12.5yum 
image show where the photometry of the peaks was obtained in all the images. 

Fig. 5.- 

Spectral energy distributions of the separate peaks in Markarian 273 as determined 
from the images presented in Figure 4. The total flux (also presented in Figure 3b) was 
measured in a 4" diameter beam centered between the peaks and was distributed between 
the two sources by doing photometry in non-overlapping boxes (displayed in the 12.5/xm 
panel of Figure 4) centered on the two peaks. The uncertainties shown are statistical. The 
systematic uncertainty depends on the wavelength, but is probably in the range 15 - 20 %. 
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Table 1 



Table 1: Basic Properties of ULIRGS Observed 



Name 


z 


log L 


Spectrum 


linear scale 










Kpc/" 


IRAS05 189-2524 


0.043 


12.08 


Sey 1 


830 


IRAS08572+3915 


0.058 


12.10 


LINER 


1200 


UGC 5101 


0.040 


12.00 


LINER 


820 


Markarian 231 


0.040 


12.52 


Sey 1 


800 


Markarian 273 


0.037 


12.13 


Sey 


800 


Arp 220 


0.018 


12.17 


LINER 


360 


IRAS17208-0011 


0.010 


12.16 


HII 


800 
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Table 2 



Table 2: Log of Keck II mid-infrared Observations of ULIRGS 



Object 


7.9fxui 




9.7fj,m 


10.3/im 


11.7//m 


12.5/xm 


17.9//m 


24.5//m 


IRAS05189 












V 




V 


IRAS08572 


V 


V 


V 






V* 


V 




UGC 5101 












V* 


V 




Mrk 231 


V 


V 


V 


V 


V 


V* 


V 




Mrk 273 


V 


V 


V 


V 


V 


V 


V 




Arp 220 


V 


V 


V 


V 


V 


V 


V 


V 


IRAS17208 


V 










V 







*size determination made with the Long Wavelength Spectrograph in March 99 
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Table 3 



Table 3: MIRLIN and IRAS Flux Densities for ULIRGS 



Object 


7.9/im 


8.8/im 


9.7/im 


10.3/im 11.7/im 
[mJy] 


12.5/im 


12.5/im* 


17.9/im 


24.5//m 


25Mm* 

[Jy] 


60//m* 


IRAS05189 










720 


740 




3.2±0.2 


3.5 


13.9 


IRAS08572 


615±18 


434±10 


221±10 




218±10 


340 


552±46 




1.8 


7.5 


UGC 5101 










185±10 


260 


256±42 




1.1 


12.1 


Mrk 231 


1370±43 


1488±18 


1081±30 


1075±21 1355±15 


1836±20 


1930 


3249±49 




8.8 


33.6 


Mrk 273 


267±15 


141±10 


92±10 


38±6 150±5 


210 ±10 


230 


442±31 




2.3 


23.7 


Arp 220 


521±15 


155±6 


124±14 


59±10 175±8 


404 ±12 


500 


1170±50 


9.8±0.2 


8.0 


104.0 


IRAS17208 


315±25 








202±13 


200 






1.7 


34.1 



The MIRLIN photometry refers to a 4" diameter beam. The uncertainties Usted in the table are statistical 
only. Photometric uncertainties in the MIRLIN data are ±5% for A < 20/im, and ±10% for A > 20/im. The 
uncertainties in the IRAS data are 5-10%. 
*Flux densities from IRAS Catalog and are not color corrected 
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Table 4 



Table 4: Properties of Compact Sources in ULIRGS 



Object 


12.5/im 


17.9;um 


24.5/im 


8.4 Ghz 


2.2 fim 




% flux+ 


size 


% flux+ 


size 


% flux+ 


size 


size 


size 




in core 


II 


in core 


II 


in core 


II 


II 


II 


IRAS05189 


100 


< 0.20 






100 


<0.26 


0.20 X 0.17 


<0.2 


IRAS08572 


~80 


<0.22* 


~100 


<0.25 






0.09 X 0.07 


<0.2 


UGC 5101 


50 


<0.25* 


60 


<0.25 






0.14 X 0.11 


<0.2 


Mrk 231 


100 


<0.13* 


85 


<0.25 






< 0.07 X 0.06 


<0.2 


Mrk 273N 


67 


<0.3 










0.32 X 0.18 


~0.3 


Mrk 273SW 


50 


<0.3 










0.18 X 0.11 


<0.2 


Arp 220 E 










25 




0.32 X 0.19 


<0.2 


Arp 220 W 










75 


0.40 


0.21 X 0.14 


0.4 


IRAS17208 


100 


2.0 










0.26 X 0.32 ^ 


2 



^percentage of flux measured in 4" diameter beam 

*size measurement from LWS observations 

^size measurement at 1.6 GHz from Martin et al. 1989 
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Table 5 



Table 5: Infrared Sizes of ULRIGs 



Object 


size of mid-ir core 


L/A[LoKpc-2] 






pc 


max 


min 


IRAS05189 


<0.2 


<166 


5.6 X 10^3 


2.1 X 10^3 


IRAS08572 


<0.22 


<260 


> 2.8 X 10^3 


> 6.0 X 10^2 


UGC 5101 


<0.25 


<205 


> 1.1 X 10^3 


> 1.0 X 10^2 


Mrk 231 


<0.13 


<104 


3.9 X 10^^ 


1.0 X 10^^ 


Mrk 273N 


<0.3 


<240 


> 2.2 X 10^3 


> 1.0 X 10^2 


Arp 220 


0.40 


140 


6.0 X 10^3 


1.2 X 10^3 


IRAS 17208 


2.0 


1600 


1.2 X 10^2 


3.0 X 10^° 
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Table 6 



Table 6: Summary of Properties of ULIRGS 



Object 


Compact core? 




Spectral Class 






MIR/NIR/Radio 


Optical 


Near IR 


ISO 


ASCA 


IRAS05189 


Y/Y/N 


Sey2 


Sey 




AGN 


IRAS08572 


Y/Y/N 


LINER 


HII 




SB 


UGC 5101 


Y/Y/Y 


LINER 




SB 


SB 


Mrk 231 


Y/Y/Y 


Seyl 


Seyl 


AGN 


SB? 


Mrk 273 


Y/Y/Y 


Sey2 


Sey? 


AGN 


AGN 


Arp 220 


Y/N/Y 


LINER 




SB 


SB? 


IRAS 17208 


N/N/? 


HII 


HII 


SB 
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